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1 Introduction 

Random walks provide a simple conventional model 
to describe various transport processes, for example 
propagation of heat or diffusion of matter through a 
medium (for a general reference see, e.g., Hughes 

(1995) ). However, in many practical cases the medium 
where the system evolves is highly irregular, due to 
factors such as defects, impurities, fluctuations etc. It 
is natural to model such irregularities as random en- 
vironment, treating the observable sample as a statis- 
tical realization of an ensemble, obtained by choos- 
ing the local characteristics of the motion (e.g., trans- 
port coefficients and driving fields) at random, ac- 
cording to a certain probability distribution. 

In the random walks context, such models are re- 
ferred to as Random Walks in Random Environments 
(RWRE). This is a relatively new chapter in applied 
probability and physics of disordered systems initi- 
ated in the 1970s. Early interest in RWRE models 
was motivated by some problems in biology, crys- 
tallography and metal physics, but later apphcations 
have spread through numerous areas (see review pa- 
pers by Alexander etal. (1981), Bouchaud and Georg( 
(1990), and a comprehensive monograph by Hughes 

(1996) ). After 30 years of extensive work, RWRE re- 
main a very active area of research, which has been 
a rich source of hard and challenging questions and 
has already led to many surprising discoveries, such 
as subdiffusive behavior, trapping effects, locaUza- 
tion, etc. It is fair to say that the RWRE paradigm 
has become firmly established in physics of random 
media, and its models, ideas, methods, results, and 
general effects have become an indispensable part of 
the standard tool kit of a mathematical physicist. 

One of the central problems in random media the- 
ory is to establish conditions ensuring homogeniza- 
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tion, whereby a given stochastic system evolving in 
a random medium can be adequately described, on 
some spatial-temporal scale, using a suitable effec- 
tive system in a homogeneous (non-random) medium. 
In particular, such systems would exhibit classical 
diffusive behavior with effective drift and diffusion 
coefficient. Such an approximation, called effective 
medium approximation (EMA), may be expected to 
be successful for systems exposed to a relatively small 
disorder of the environment. However, in certain cir- 
cumstances EMA may fail due to atypical environ- 
ment configurations ("large deviations") leading to 
various anomalous effects. For instance, with small 
but positive probability a reaUzation of the environ- 
ment may create "traps" that would hold the particle 
for anomalously long time, resulting in the subdif- 
fusive behavior, with the mean square displacement 
growing slower than linearly in time. 

RWRE models have been studied by various non- 
rigorous methods including Monte Carlo simulations, 
series expansions, and the renormalization group tech- 
niques (see more details in the above references), but 
only few models have been analyzed rigorously, es- 
pecially in dimensions greater than one. The situa- 
tion is much more satisfactory in the one-dimensional 
case, where the mathematical theory has matured and 
the RWRE dynamics has been understood fairly well. 

The goal of this article is to give a brief introduc- 
tion to the beautiful area of RWRE. The principal 
model to be discussed is a random walk with nearest- 
neighbor jumps in independent identically distributed 
(i.i.d.) random environment in one dimension, al- 
though we shall also comment on some generaliza- 
tions. The focus is on rigorous results; however, 
heuristics will be used freely to motivate the ideas 
and explain the approaches and proofs. In a few 
cases, sketches of the proofs have been included, which 
should help appreciate the flavor of the results and 
methods. 

1.1 Ordinary Random Walks: A Reminder 

To put our exposition in perspective, let us give a 
brief account of a few basic concepts and facts for 
ordinary random walks, that is, evolving in a non- 
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random environment (see further details in Hughes 
1995). In such models, space is modelled using a 
suitable graph, e.g., a d-dimensional integer lattice 
U^, while time may be discrete or continuous. The 
latter distinction is not essential, and in this article we 
will mostly focus on the discrete-time case. The ran- 
dom mechanism of spatial motion is then determined 
by the given transition probabilities (probabilities of 
jumps) at each site of the graph. In the lattice case, it 
is usually assumed that the walk is translation invari- 
ant, so that at each step distribution of jumps is the 
same, with no regard to the current location of the 
walk. 

In one dimension (d = 1), the simple (nearest- 
neighbor) random walk may move one step to the 
right or to the left at a time, with some probabilities 
p and q = 1 — p, respectively. An important as- 
sumption is that only the current location of the walk 
determines the random motion mechanism, whereas 
the past history is not relevant. In terms of proba- 
bility theory, such a process is referred to as Markov 
chain. Thus, assuming that the walk starts at the ori- 
gin, its position after n steps can be represented as 
the sum of consecutive displacements, X„ = Zi + 
• • • + Zn, where Zi are independent random vari- 
ables with the same distribution P{Zi = 1} = p, 
P{Zi = -1} = q. 

The strong law of large numbers (LLN) states that 
almost surely (i.e., with probability 1) 

lim — = EZi=p — P-a.s. (1) 

n— »oo n 

where E denotes expectation (mean value) with re- 
spect to P. This result shows that the random walk 
moves with the asymptotic average velocity close to 
p — q. It follows that if p — q ^ then the pro- 
cess Xn, with probability 1, will ultimately drift to 
infinity (more precisely, +oo if p — q > and — oo 
if p — q < 0). In particular, in this case the ran- 
dom walk may return to the origin (and in fact visit 
any site on Z) only finitely many times. Such be- 
havior is called transient. However, in the symmet- 
ric case (i.e., p = q = 0.5) the average velocity 
vanishes, so the above argument fails. In this case 
the walk behavior appears to be more complicated, 
as it makes increasingly large excursions both to the 
right and to the left, so that lim„^oo = +oo, 
lim^^^ Xn = — oo (P-a.s.). This implies that a 
symmetric random walk in one dimension is recur- 
rent, in that it visits the origin (and indeed any site 



on Z) infinitely often. Moreover, it can be shown 
to be null-recurrent, which means that the expected 
time to return to the origin is infinite. That is to say, 
return to the origin is guaranteed, but it takes very 
long until this happens. 

Fluctuations of the random walk can be charac- 
terized further via the central limit theorem (CLT), 
which amounts to saying that the distribution of X„ 
is asymptotically normal, with mean n{p — q) and 
variance Anpq: 

n^oo [ ^J'inpq J 

= ^(x):=^ r e-y'l'^dy. (2) 

V ^TT J-oo 

These results can be extended to more general walks 
in one dimension, and also to higher dimensions. For 
instance, the criterion of recurrence for a general one- 
dimensional random walk is that it is unbiased, E{Xi~ 
Xq) = 0. In the two-dimensional case, in addition 
one needs E\Xi — Xq\'^ < oo. In higher dimen- 
sions, any random walk (which does not reduce to 
lower dimension) is transient. 

1.2 Random Environments and Random 
Walks 

The definition of an RWRE involves two ingredi- 
ents: (i) the environment, which is randomly cho- 
sen but remains fixed throughout the time evolution, 
and (ii) the random walk, whose transition proba- 
bilities are determined by the environment. The set 
of environments (sample space) is denoted by $7 = 
{uj}, and we use P to denote the probability distri- 
bution on this space. For each u; G Q, we define 
the random walk in the enviroimient u as the (time- 
homogeneous) Markov chain {Xt, t = 0, 1, 2, ... } 
on Z*^ with certain (random) transition probabiMties 

p{x,y,uj) = P^{X, = y\Xo = x}. (3) 

The probability measure P"^ that determines the dis- 
tribution of the random walk in a given environment 
Lo is referred to as the quenched law. We often use a 
subindex to indicate the initial position of the walk, 
so that e.g. P%{Xq = x} = 1. 

By averaging the quenched probability fur- 
ther, with respect to the environment distribution, we 
obtain the annealed measure Px = P x P^, which 
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determines the probability law of the RWRE: 

P,iA)= [ P-(A)P(da;) = EP-(^). (4) 
Jn 

Expectation with respect to the annealed measure Px 
will be denoted by E^^^. 

Equation (4) implies that if some property A of the 
RWRE holds almost surely (a.s.) with respect to the 
quenched law for almost all environments (i.e., 
for all Lu e ^' such that F{^}') = 1), then this prop- 
erty is also true with probabiUty 1 under the annealed 
law Px- 

Note that the random walk Xn is a Markov chain 
only conditionally on the fixed envirormient (i.e., with 
respect to P^), but the Markov property fails under 
the annealed measure Px- This is because the past 
history caimot be neglected, as it tells what informa- 
tion about the medium must be taken into account 
when averaging with respect to environment. That 
is to say, the walk learns more about the environ- 
ment by taking more steps. (This idea motivates the 
method of "environment viewed from the particle", 
see Section 7 below.) 

The simplest model is the nearest-neighbor one- 
dimensional walk, with transition probabilities 

' Px if y = x-l- 1, 
p{x,y,u) = < Qx ify = x-l, 
otherwise, 

where px and qx = 1 — Px (x ^ li) are random vari- 
ables on the probability space ($7, P). That is to say, 
given the environment a; G the random walk cur- 
rently at point a; e Z will make a one-unit step to the 
right, with probability px, or to the left, with prob- 
ability Qx- Here the environment is determined by 
the sequence of random variables {px}- For the most 
of the article, we assume that the random probabil- 
ities {px, X G Z} are independent and identically 
distributed (i.i.d.), which is referred to as i.i.d. en- 
vironment. Some extensions to more general envi- 
ronments will be mentioned briefly in Section 9. The 
study of RWRE is simplified under the following nat- 
ural condition called (uniform) ellipticity: 

0<(5<p^<l-(5<l, xGZ, P-a.s. (5) 

which will be frequently assumed in the sequel. 

2 Transience and Recurrence 

In this section, we discuss a criterion for the RWRE 
to be transient or recurrent. The following theorem 



is due to Solomon (1975). 

Theorem 1. Set px ■= qx/Px> x G Z, and r) := 

Eln/90- 

(i) If T] ^ then Xt is transient (Po-a.s.); more- 
over, if r] < then limt_>oo Xf = -|-oo, while if 
r] > then limt^oo Xt = — oo (Po-a.s.). 

(ii) If 7] = then Xt is recurrent (Po-a.s.); more- 
over, 

lim Xt = +00, lim Xt = —00, Po-a.s. 

Let us sketch the proof. Consider the hitting times 
Tx := min{t > : Xt = x} and denote by fxy the 
quenched first-passage probability from x to y: 

fxy ■= P%{1 <Ty< 00}. 

Starting from the first step of the walk may be ei- 
ther to the right or to the left, hence by the Markov 
property the return probability /oo can be decom- 
posed as 

/oo = Po/10 + qof-1,0- (6) 
To evaluate /lo, forn > 1 set 

Ui = «f ^ := Pt{To <Tx}, 0<i<x, 

which is the probability to reach prior to x, starting 
from i. Clearly, 



/lo = lim u^^K 

x—*oo 



(7) 



Decomposition with respect to the first step yields 
the difference equation 

Ui = piUi+i -\- qiUi-i, <i <n, (8) 

with the boundary conditions 

no = 1, ux = 0. (9) 

Using Px + Qx = I5 eqn (8) can be rewritten as 

Ux+l -Ux = PxiUx - Ux-l), 

whence by iterations 

x 

Ux+l -Ux = (ui - no) W Pj. (10) 

i=i 

Summing over x and using the boundary conditions 
(9) we obtain 

/n— 1 X \ 



Va;=0 j=l 
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(if X = 0, the product over j is interpreted as 1). In 
view of eqn (7) it follows that /lo = 1 if and only if 
the right-hand side of eqn (11) tends to 0, that is, 

OO X 

J2 exp(y^) = OO, := ^ Inp^-. (12) 

x=l j=l 

Note that the random variables In pj are i.i.d., hence 
by the strong LLN 



a 



lim — = E In po 

a;— >oo X 



-a.s. 



That is, the general term of the series (12) for large x 
behaves like exp(x?7), hence for 77 > the condition 
(12) holds true (and so /lo = 1), whereas for r/ < 
it fails (and so /lo < 1). 

By interchanging the roles of and q^, we also 
have /-i,o < 1 if ?? > and /-i,o = 1 if r/ < 0. 
From eqn (6) it then follows that in both cases /oo < 
1, i.e. the random walk is transient. 

In the critical case, r/ = 0, by a general result from 
probability theory, > for infinitely many x (P- 
a.s.), and so the series in eqn (12) diverges. Hence, 
/lo = 1 and, similarly, /_i,o = 1, so by eqn (6) 
/oo = 1, i.e. the random walk is recurrent. 

It may be surprising that the critical parameter ap- 
pears in the form 77 = E In po, as it is probably more 
natural to expect, by analogy with the ordinary ran- 
dom walk, that the RWRE criterion would be based 
on the mean drift, E(po — <?o)- In the next section we 
will see that the sign of d may be misleading. 

A canonical model of RWRE is specified by the 
assumption that the random variables px take only 
two values, j3 and I — (3, with probabilities 

= P} = a, Pfe = 1 - /?} = 1 - a, (13) 

where 0<a<l, 0</3<l. Here 77 = (2a - 
l)ln(l -I- (1 - 2/3) /P), and it is easy to see that, 
e.g., r/ < if a < ^, P < ^ or a > ^, P > 5. 
The recurrent region where rj = sphts into two 
lines, /5 = 5 and a = ^. Note that the first case 
is degenerate and amounts to the ordinary symmetric 
random walk, while the second one (except where 
P = ^) corresponds to Sinai's problem (see Section 
6). A "phase diagram" for this model, showing vari- 
ous limiting regimes as a function of the parameters 
a, P, is presented in Figure 1 . 
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Figure 1 : Phase diagram for the canonical model, eqn 
(13). In the regions where < or 77 > 0, the RWRE is 
transient to -|-oo or —00, respectively. The recurrent case, 
r] = 0, arises when a = i or /? = ^ . The asymptotic ve- 
locity V := lim(_^oc Xi /tAs given by eqn (14). Adapted 
from Hughes B.D. (1996) Random Walks and Random 
Environments. Volume 2: Random Environments, Ch. 6, 
p. 39 1 . Clarendon Press, Oxford, by permission of Oxford 
University Press. 

3 Asymptotic Velocity 

In the transient case the walk escapes to infinity, and 
it is reasonable to ask at what speed. For a non- 
random environment, px = p, the answer is given by 
the LLN, eqn (1). For the simple RWRE, the asymp- 
totic velocity was obtained by Solomon (1975). Note 
that by Jensen's inequahty, (E po)~^ < E p^^. 



Theorem 2. The limit v 
( Po-a.s.) and is given by 

1-Epo 
1 + Epo 



limt_>oo^t/i exists 



V = i 



ifEpo < 1, 



otherwise. 



(14) 



Thus, the RWRE has a well-defined non-zero 

asymptotic velocity except when {'Epo)~^ < 1 < 
E/9q^. For instance, in the canonical example eqn 
(13) (see Figure 1) the criterion E po < 1 for the ve- 
locity V to be positive amounts to the condition that 
both (1 — a)/a and (1 — P)/p\ie on the same side 
of point 1. 

The key idea of the proof is to analyze the hitting 
times Tn first, deducing results for the walk Xt later. 
More specifically, set n = Ti — Ti-i, which is the 
time to hit i after hitting i—1 (providing that i > Xq). 
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If Xo = and n > 1 then T„ = n + • • • + r„. 
Note that in fixed environment uj the random vari- 
ables {ti} are independent, since the quenched ran- 
dom walk "forgets" its past. Although there is no 
independence with respect to the annealed probabil- 
ity measure Pq, one can show that, due to the i.i.d. 
property of the environment, the sequence {rj} is er- 
godic and therefore satisfies the LLN: 



n 



n 



In turn, this implies 

Xt ^ 
t ^ 



Eon' 



Eon, Po-a.s. 



Po-a.s. (15) 



(the clue is to note that Xt„ = n). 

To compute the mean value Eq t\, observe that 



Furthermore, by Jensen's inequality 

E/9o = Epo'-l > (Epo)~'-l, 
so eqn (14) implies that if E po < 1 then 

0<i;<2E;jo-l = E(po-go), 

and the inequality is strict if po is genuinely ran- 
dom (i.e., does not reduce to a constant). Hence, 
the asymptotic velocity v is less than the mean drift 
E(po — qo), which is yet another evidence of slow- 
down. What is even more surprising is that it is pos- 
sible to have E(po — 9o) > but 77 = Elnpo > 0, 
so that Po-a.s. Xt —00 (although with velocity 
V = 0). Indeed, following Sznitman (2004) suppose 
that 



7} = 1 - a, 



L{Xi=l} 



+ 1 



{^1 



i}(1 + ^o + tO 



(16) 



where 1a is the indicator of event A and Tq, t[ are, 
respectively, the times to get from —1 to and then 
from to 1. Taking expectations in a fixed environ- 
ment UJ, we obtain 



with a> \. Then Epo > a/3 > | if 1 > /9 > ^, 
hence E(po — <?o) = 2Epo — 1 > 0. On the other 
hand, 

1-/3 1 - ^ 
Eln/9o = Q;ln— ^-F(l-a)ln > 0, 

P 7 



n = Po + go ( 1 + £^0 ^0 + ^0 n ) , (17) if 7 is sufficiently small. 



and so 



^^n = l + po + Po^o^, 



0- 



(18) 



Note that Eq Tq is a function of {p^ , x < 0} and 
hence is independent of po = qo/po- Averaging eqn 
(18) over the environment and using Eq Tq = Eq t\ 
yields 



EqTx 



cxD if Epo > Ij 



and by eqn (15) "half" of eqn (14) follows. The other 
half, in terms of E/3q^, can be obtained by inter- 
changing the roles of px and Qx, whereby po is re- 
placed with Pq^. 

Let us make a few remarks concerning Theorems 
1 and 2. First of all, note that by Jensen's inequality 
E Inpo < InE po5 with a strict inequality whenever 
Po is non-degenerate. Therefore, it may be possible 
that, with Po-probability 1, ^ 00 but X^/t ^ 
(see Figure 1). This is quite unusual as compared to 
the ordinary random walk (see Section 1.1), and in- 
dicates some kind of slowdown in the transient case. 



4 Critical Exponent, Excursions and 
Traps 

Extending the previous analysis of the hitting times, 
one can obtain useful information about the limit dis- 
tribution of Tn (and hence Xt). To appreciate this, 
note that from the recursion (16) it follows 

rf = l{Xi=i} + l{Xi=-i}(l + ro + rO", 

and, similarly to eqn (17), 

E^Tl=po + qoE'^{l + T'o + T[y. 

Taking here expectation E, one can deduce that Eq rf 
< 00 if and only if E Pq < 1. Therefore, it is natural 
to expect that the root k of the equation 



Ep« = l 



(20) 



plays the role of a critical exponent responsible for 
the growth rate (and hence, for the type of the limit 
distribution) of the sum T„ = n H hr^. In partic- 
ular, by analogy with sums of i.i.d. random variables 
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one can expect that if k > 2 then T„ is asymptot- 
ically normal, with the standard scaling ^Jn, while 
for K < 2 the limit law of Tn is stable (with index k) 
under scaling f» n^/*^. 

Alternatively, eqn (20) can be obtained from con- 
sideration of excursions of the random walk. Let T^y 
be the left excursion time from site 1, that is the time 
to return to 1 after moving to the left at the first step. 
If ?7 = Elnpo < 0, then r/i[ < cx) (Po-a.s.). Fixing 
an environment w, let u^i = T/^ be the quenched 
mean duration of the excursion and observe that 
w\ = \ ^ T\, where t\ is the time to get back to 
1 after stepping to 0. 

As a matter of fact, this representation and eqn 
(19) imply that the annealed mean duration of the 
left excursion, T^x, is given by 

TO J \ \ ifEpo<l, ,r,.. 

Eu;i = < 1 — Epo (21) 
\ oo ifEpo>l- 

Note that in the latter case (and bearing in mind r] < 
0), the random walk starting from 1 will eventually 
drift to -l-oo, thus making only a finite number of 
visits to 0, but the expected number of such visits is 
infinite. 

In fact, our goal here is to characterize the distri- 
bution of wi under the law P. To this end, observe 
that the excursion T^^ involves at least two steps (the 
first and the last ones) and, possibly, several left ex- 
cursions from 0, each with mean time u^o = Eq Tq^. 
Therefore, 

oo 

u;i = 2 -I- ^ qipoijwo) = 2 + pqwq. (22) 

i=i 

By the translation invariance of the environment, the 
random variables wi and wq have the same distri- 
bution. Furthermore, similarly to recursion (22), we 
have Wq = 2 + p-iw-i. This implies that wq is 
a function of Px with x < — 1 only, and hence u^o 
and pq are independent random variables. Introduc- 
ing the Laplace transform (j){s) = E exp(— swi) and 
conditioning on pQ, from eqn (22) we get the equa- 
tion 



E0(spo)- 



(23) 



Suppose that 



1 - (/.(s) 



~ as 



Expanding the product on the right, one can see that 
a solution with k = 1 is possible only if E po < in 
which case 

9 



then eqn (23) amounts to 



We have already obtained this result in eqn (21). 

The case k < 1 is possible if E Pq = 1, which 
is exactly eqn (20). Returning to wi, one expects a 
slow decay of the distribution tail, 

P{u;i > ~ t^oo. 

In particular, in this case the aimealed mean duration 
of the left excursion appears to be infinite. 

Although the above considerations point to the crit- 
ical parameter n, eqn (20), which may be expected 
to determine the slowdown scale, they provide little 
explanation of a mechanism of the slowdown phe- 
nomenon. Heuristically, it is natural to attribute the 
slowdown effects to the presence of traps in the envi- 
ronment, which may be thought of as regions that are 
easy to enter but hard to leave. In the one-dimensional 
case, such a trap would occur, for example, between 
two long series of successive sites where the proba- 
biUties px are fairly large (on the left) and small (on 
the right). 

Remarkably, traps can be characterized quantita- 
tively with regard to the properties of the random en- 
viroimient, by linking them to certain large deviation 
effects (see Sznitman (2002, 2004)). The key role 
in this analysis is played by the function F{u) := 
InEp^, M G M. Suppose that = Elnpo < (so 
that by Theorem 1 the RWRE tends to -|-co, Po-a.s.) 
and also that E po > 1 and E pQ^ > 1 (so that by 
Theorem 2, v = 0). The latter means that F{1) > 
and F{—1) > 0, and since F is a smooth strictly 
convex function and F(0) = 0, it follows that there 
is the second root < k < 1, so that F{k) = 0, i.e., 
E p(j = 1 (cf. eqn (20)). 

Let us estimate the probability to have a ttap in 
U = [—L,L] where the RWRE will spend anoma- 
lously long time. Using eqn (11), observe that 

Pfin <Tl+i}>1- exp{-LSL}, 

where Sl ■= J2x=i^''^ Px i] < as L —>^ 
oo. However, due to large deviations Sl may exceed 
level e > with probability 



l-as'' + --- = {l-2s + ---){l-as^Ep!^ + ---). P{5l > e} ~ exp{-L/(e)}, L 



oo. 
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where I{x) := sup^jux — F{u)} is the Legendre 
transform of F. We can optimize this estimate by 
assuming that eL > Inn and minimizing the ratio 
/(e) /e. Note that F{u) can be expressed via the 
inverse Legendre transform, F{u) = s\ip^{xu — 
I{x)}, and it is easy to see that if k := mine>o /(e) / e 
then F{k) = 0, so k is the second (positive) root 
of F. 

The "left" probability P'^^JTq < T^l-i} is esti- 
mated in a similar fashion, and one can deduce that 
for some constants K > 0, c > and any k' > k, 
for large n 

p|p^|max|Xfc| < Klnnj > c} > n'^' . 

That is to say, this is a bound on the probability to 
see a trap centered at 0, of Size ^ Inn, which will 
retain the RWRE for at least time n. It can be shown 
that, typically, there will be many such traps both 
in [— n*^ , 0] and [0, n'^ ], which will essentially pre- 
vent the RWRE from moving at distance n'^ from 
the origin before time n. In particular, it follows that 
lim^^oo Xn/n'^ = for any k' > k, so recalling 
that < K < 1, we have indeed a sublinear growth 
of Xn. This result is more informative as compared 
to Theorem 2 (the case v = 0), and it clarifies the 
role of traps (see more details in Sznitman (2004)). 
The non-trivial behavior of the RWRE on the precise 
growth scale, n**, is characterized in the next section. 

5 Limit Distributions 

Considerations in Section 4 suggest that the exponent 

K, defined as the solution of eqn (20), characterizes 
environments in terms of duration of left excursions. 
These heuristic arguments are confirmed by a limit 
theorem by Kesten et al. (1975), which specifies the 
slowdown scale. We state here the most striking part 
of their result. Denote In^n := maxjln n, 0}; by an 
arithmetic distribution one means a probability law 
on M concentrated on the set of points of the form 0, 
±c, ±2c, . . . 

Theorem 3. Assume that — oo<r7 = Eln/9o<0 
and the distribution of In po is non-arithmetic ( ex- 
cluding a possible atom at — ocj. Suppose that the 
root K of equation (20) is such that < k < 1 and 



E Pq In^po < oo. Then 

lim Po{n-^'''Tn<t} = LS), 

n— »oo 

lim Po{r''Xt <x} = l- L4x-V«), 

t— >oo 

where //«(•) is the distribution function of a stable 
law with index k, concentrated on [0, oo). 

General information on stable laws can be found 
in many probability books; we only mention here 
that the Laplace transform of a stable distribution on 
[0, oo) with index k has the form 0(s) = exp{— Cs''} 

Kesten et al. (1975) also consider the case k > 1. 
Note that for k > 1, we have E po < {E Pq)^/'^ = 1, 
so u > by eqn (14). For example, if k > 2 then, as 
expected (see Section 4), 

lim pJ^:^<t\ = m, 

Let us describe an elegant idea of the proof based 
on a suitable renewal structure, (i) Let ?7" (i < n) 
be the number of left excursions starting from i up to 
time Tn, and note that r„ = n + 2 ^ • U'^. Since the 
walk is transient to +oo, the sum X^j<o f^f is finite 
(Po-a.s.) and so does not affect the limit, (ii) Ob- 
serve that if the environment lv is fixed then the con- 
ditional distribution of ?7", given • • • ' — 0' 
is the same as the distribution of the sum of 1 + 
J7JYi i-i-d. random variables Vi , V2 , . . . , each with 
geometric distribution -Pgi^ = k} = pjqj (k = 
0,1,2,...). Therefore, the sum Y17=i (read from 
right to left) can be represented as Ylt=o where 
Zq = 0, Zi, Z2, . . . is a branching process (in ran- 
dom environment {pj}) with one immigrant at each 
step and the geometric offspring distribution with pa- 
rameter Pj for each particle present at time j. (iii) 
Consider the successive "regeneration" times t^, at 
which the process Zt vanishes. The partial sums 
Wk := X]r*<t<r* form an i.i.d. sequence, and 
the proof amounts to showing that the sum of Wk 
has a stable limit of index k. (iv) Finally, the dis- 
tribution of Wo can be approximated using Mq := 
X^t^i nj=o Pj (11))' which is the quenched 

mean number of total progeny of the immigrant at 
time t = 0. Using Kesten's renewal theorem, it can 
be checked that P{Mo > x} ~ Kx"'^ as x ^ 00, 
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so Mq is in the domain of attraction of a stable law 
with index k, and the result follows. 

Let us emphasize the significance of the regenera- 
tion times T* . Returning to the original random walk, 
one can see that these are times at which the RWRE 
hits a new "record" on its way to +00, never to back- 
track again. The same idea plays a crucial role in 
the analysis of the RWRE in higher dimensions (see 
Sections 10.1, 10.2 below). 

Finally, note that the condition —00 < 77 < al- 
lows P{po = 1} > 0, so the distribution of po may 
have an atom at (and hence In po at —00). In view 
of eqn (20), no atom is possible at +00. The restric- 
tion for the distribution of In po to be non-arithmetic 
is important. This will be illustrated in Section 8 
where we discuss the model of random diodes. 

6 Sinai's Localization 

The results discussed in Section 5 indicate that the 
less transient the RWRE is (i.e., the critical exponent 
decreasing to zero), the slower it moves. Sinai (1982) 
proved a remarkable theorem showing that for the 
recurrent RWRE (i.e., with = E In po = 0), the 
slowdown effect is exhibited in a striking way. 

Theorem 4. Suppose that the environment {px\ is 
i.i.d. and elliptic, eqn (5), and assume that E,ln po = 
0, with P{po = !}<!. Denote := ElnVo, < 
(T^ < 00. Then there exists a function Wn = Wn{oj) 
of the random environment such that for any e > 



lim Po 

n— >oo 



\r? n 



> e 



0. 



Moreover, Wn has a limit distribution: 
lim P{Wn <x} = G{x), 



(24) 



(25) 



and thus also the distribution ofa^Xn/ In^ n under 
Pq converges to the same distribution G{x). 

Sinai's theorem shows that in the recurrent case, 
the RWRE considered on the spatial scale In^ n be- 
comes localized near some random point (depending 
on the environment only). This phenomenon, fre- 
quently referred to as Sinai's localization, indicates 
an extremely strong slowdown of the motion as com- 
pared with the ordinary diffusive behavior. 

Following Revesz (1990), let us explain heuristi- 



eqn (11) as 

/ n-l . -1 

P^iTn < To} = (1 + ''MYx)j , (26) 

where is defined in (12). By the central limit the- 
orem, the typical size of \Yx\ for large x is of order 
of ^/x, and so eqn (26) yields 

P'i{Tn<To}^exp{-^/^}. 

This suggests that the walk started at site 1 will make 
about ex^p{^/n} visits to the origin before reaching 
level n. Therefore, the first passage to site n takes 
at least time ^ exp{i/n }. In other words, one may 
expect that a typical displacement after n steps will 
be of order of In^ n (cf. eqn (24)). This argument 
also indicates, in the spirit of the trapping mechanism 
of slowdown discussed at the end of Section 4, that 
there is typically a trap of size In^ n, which retains 
the RWRE until time n. 

It has been shown (independently by H. Kesten 
and A.O. Golosov) that the Umit in (25) coincides 
with the distribution of a certain functional of the 
standard Brownian motion, with the density function 



G'(x) = -yif^exp 
k=0 



(2A; + 1)V 



7 Environment Viewed from the 
Particle 

This important technique, dating back to Kozlov and 
Molchanov (1984), has proved to be quite efficient in 
the study of random motions in random media. The 
basic idea is to focus on the evolution of the environ- 
ment viewed from the current position of the walk. 

Let 9 be the shift operator acting on the space of 
environments $7 = {uo} as follows: 



^ = {Px} 1-^0) = {Px-l}- 

Consider the process 

which describes the state of the environment from 
the point of view of an observer moving along with 
the random walk X„. One can show that u;„ is a 
Markov chain (with respect to both Pq and Pq), with 
the transition kernel 



cally why X„ is measured on the scale In n. Rewrite T{uj, dco') = po dgoj{d^') + qo ^e-^A^^') (27) 
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and the respective initial law 5^^ or P (here 5^ is the 
Dirac measure, i.e., unit mass at w). 

This fact as it stands may not seem to be of any 
practical use, since the state space of this Markov 
chain is very complex. However, the great advantage 
is that one can find an explicit invariant probability 
Q for the kernel T (i.e., such that Q T = Q), which 
is absolutely continuous with respect to P. 

More specifically, assume that E po < 1 set 
Q = j\uj) P, where (cf. eqn (14)) 

1-Epo 



Due to the Markov property, the process M„ := X„— 
Dn is a martingale with respect to the natural filtra- 
tion Tn = (^{Xi, . . . , Xn} and the law Pq, 

E'^[Mn+l\Tn]= Mn (P^a.S.), 

and it has bounded jumps, |M„ — -M„_i| < 2. By 
general results, this implies Mn/n —>■ {Pq-sl.s.). 
On the other hand, by Birkhoff 's ergodic theorem 



lim 

n—*oo n 



d(0,w)Q(du;), Po-a.s. 



f = v{l + po) ^ n Pj' ^ = ^ — E^' ^^^^ ^^^^ integral is easily evaluated to yield 



x=0 j=l 

Using independence of {px}, we note 



„ oo 

/ Q(du;)=E/ = (l-Epo) V(Epor = 1, 

hence Q is a probability measure on O. Furthermore, 
for any bounded measurable function 3 on O we have 

qTg= f Tg{Lo)Q{du;)=EfTg 
Jn 

= E{f[poigoe) + qo{goe-')]} ^^^^ 
= E{g[{pof)oe-^ + {qof)o9)]. 

By eqn (28), 

00 X 
x=Q j=l 

(00 X \ 
x=Oi=i ^ 



and similarly 



{qof) oe = -v + 



1 + Po 

So from eqn (29) we obtain 

qTg = E{gf)= [ g{co)Q{dLo) = Qg, 
Jn 

which proves the invariance of Q. 

To illustrate the environment method, let us sketch 
the proof of Solomon's result on the asymptotic ve- 
locity (see Theorem 2 in Section 3). Set d{x, to) := 
- Xo) = Px- Qx- Noting that d{x,uj) = 
d(0,6'^w), define 

n n 



Hpo -qo)f = v-EY,Yl - Po) 

x=Oj=l 



vil-Epo)J2i^Po) 



x=0 



and the first part of the formula (14) follows. 

The case Epo > 1 can be handled using a com- 
parison argument (Sznitman 2004). Observe that if 
Px < Px for all X then for the corresponding random 
walks we have Xt < Xt (-Pq "^-S-)- We now define a 
suitable dominating random medium by setting (for 
7>0) 



Px 



Px 



+ 



7 



> Px 



1+7 1+7 

Then E po = E qo/ {po + j) < 1 if 7 is large enough, 
so by the first part of the theorem, Pq-sls., 



- — Xfi 

lim < lim 

n— »oo n n— +00 n 



Epo 



(30) 



1+Epo 

Note that E po is a continuous function of 7 with val- 
ues in [0, E Po] 9 1» so there exists 7* such that E po 
attains the value 1. Passing to the limit in eqn (30) 
as 7 t 7*5 we obtain lim„^oo Xn/n < (Po'^-S-)- 
Similarly, we get the reverse inequaUty, which proves 
the second part of the theorem. 

A more prominent advantage of the environment 
method is that it naturally leads to statements of CLT 
type. A key step is to find a function H{x, t, uj) = 
X — vt + h(x, iu) (called harmonic coordinate) such 
that the process H{Xn, n, uj) is a martingale. To this 
end, by the Markov property it suffices to have 

E'^^H{Xn+i,n + l,iv) = H{Xn,n,uj), P^ -a.s. 

For A(x, uj) := h{x + 1,uj) — h{x, u) this condition 
leads to the equation 

A{x,uj) = Px^ix -l,uj) + v-l + {l + v)px. 
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If E po < 1 (so that V > 0), there exists a bounded 
solution 

oo k 
k=0 i=0 

and we note that A{x, u) = A(0, O^u) is a station- 
ary sequence with mean IEA(x,a;) = 0. Finally, 
setting h{0, w) = we find 



/■ x-l 



h{x,uj) = < 



^A(fc,cc;), x>0, 
fc=o 

— X 

^A(-fc,u;), X < 0. 



fe=i 



As a result, we have the representation 



nv 



H{Xn, n, iv) + h{Xn, LO) 



(31) 



For a fixed to, one can apply a suitable CLT for mar- 
tingale differences to the martingale term in (31), 
while using that Xn ~ nv {Po-a.s.), the second term 
in (31) is approximated by the sum Ylk=o ^i^i '^)' 
which can be handled via a CLT for stationary se- 
quences. This way, we arrive at the following result. 

Theorem 5. Suppose that the environment is elliptic, 
eqn (5), and such that E Pg^^ < 1 for some e > 
(which implies that K pq < land hence v > 0). Then 
there exists a non-random cr^ > such that 



lim Pq 

n—^oo 



Xr, 



nv 



< X 



Note that this theorem is parallel to the result by 
Kesten et al. (1975) on asymptotic normality when 
K > 2 (see Section 5). The assumptions in Theorem 
5 as stated are a bit more restrictive than in Theorem 
3, but they can be relaxed. More importantly, the 
environment method proves to be quite efficient in 
more general situations, including non-i.i.d. environ- 
ments and higher dimensions (at least in some cases, 
e.g., for random bonds RWRE and balanced RWRE. 

8 Diode Model 

In the preceding sections (except in Section 5, where 
however we were limited to a non-arithmetic case), 
we assumed that < < 1 and therefore excluded 
the situation where there are sites through which mo- 
tion is permitted in one direction only. Allowing for 



such a possibility leads to the diode model (Solomon 
1975). Specifically, suppose that 

= /?} = «, Pfe = l} = l-a, (32) 

with 0<q:<1, 0</3<1, so that with probabil- 
ity a a point x G Z is a usual two-way site and with 
probability 1 — a it is a repelling barrier ("diode"), 
through which passage is only possible from left to 
right. This is an interesting example of statistically 
inhomogeneous medium, where the particle motion 
is strongly irreversible due to the presence of special 
semi-penetrable nodes. The principal mathematical 
advantage of such a model is that the random walk 
can be decomposed into independent excursions from 
one diode to the next. 

Due to diodes the random walk will eventually 
drift to -l-oo. If /? > \, then on average it moves 
faster than in a non-random environment with px = 
j3. The situation where /? < g is potentially more in- 
teresting, as then there is a competition between the 
local drift of the walk to the left (in ordinary sites) 
and the presence of repelling diodes on its way. Note 
that E po = ctp, where p := (1 — /3)//3, so the con- 
dition E po < 1 amounts to /? > a/(l + a). In this 
case (which includes (3 > ^), formula (14) for the 
asymptotic velocity applies. 

As explained in Section 4, the quenched mean du- 
ration w of the left excursion has Laplace transform 
given by eqn (23), which now reads 

4>{s) = e"^''{l - a + a(j){sp)}. 

This equation is easily solved by iterations: 

oo 

0(s) = (l-a)^aV^*^ 



(33) 



k=0 
k 

i=o 

hence the distribution of w is given by 

F{w = tk} = {l-a)a'', A; = 0,1,... 

This result has a transparent probabilistic meaning. 

In fact, the factor ( 1 — a) a'^ is the probability that the 
nearest diode on the left of the starting point occurs 
at distance k + 1, whereas tk is the corresponding 
mean excursion time. Note that formula (33) for 
easily follows from the recursion = 2 + pt^-i (cf. 
eqn (22)) with the boundary condition to = 2. 
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A self-similar hierarchy of time scales (33) indi- 
cates that the process will exhibit temporal oscilla- 
tions. Indeed, for ap > 1 the average waiting time 
until passing through a valley of ordinary sites of 
length k is asymptotically proportional to ~ 2p^, 
so one may expect the annealed mean displacement 
£^0 to have a local minimum att ^ t^. Passing to 
logarithms, we note that In tk+i — In tk ~ In p, which 
suggests the occurrence of persistent oscillations on 
the logarithmic time scale, with period In p. This was 
confirmed by Bemasconi and Schneider (1985) who 
showed that for ap > 1 

Eq Xn ~ n''F(lnn), n^oo, (34) 

where k = — Ina/lnp < 1 is the solution of eqn 
(20) and the function F is periodic with period In p 
(see Figure 2). 

0.21 ^ 1 ' 1 1 1 1 1 1 1 '~ 

-r ■ ' • - 

T . ... ..• V 

1 I li I 1 <i ' it. I L— 

D 2 1 6 e 10 

In t 

Figure 2: Temporal oscillations for the diode model, eqn 
(32). Here a = 0.3 and p = 1/0.09, so that ap > 1 
and K = \. The dots represent an average of Monte 
Carlo simulations over 10 000 samples of the environment 
with a random walk of 200 000 steps in each realization. 
The broken curve refers to the exact asymptotic solution 
(34). The arrows indicate the simulated locations of the 
minima tk, the asymptotic spacing of which is predicted 
to be Inp « 2.41. Reproduced from Bernasconi J. and 
Schneider W.R. (1982) Diffusion on a one-dimensional 
lattice with random asymmetric transition rates. Journal 
of Physics A: Mathematical and General, 15, L729-L734, 
by permission of lOP Publishing Ltd. 

In contrast, for ap = 1 one has 



and there are no oscillations of the above kind. 

These results illuminate the earlier analysis of the 
diode model by Solomon (1975), which in the main 
has revealed the following. If ap = 1 then X„ satis- 
fies the strong LLN: 

Xn Inp 
hm — = Po-a.s., 



while in the case ap > 1 the asymptotic behavior of 
Xn is quite complicated and unusual: if ^ oo is 
a sequence of integers such that {Inni} 7 (here 
{a} = a — [a] denotes the fractional part of a), then 
the distribution of n^'^Xm under Pq converges to 
a non-degenerate distribution which depends on 7. 
Thus, the very existence of the limiting distribution 
of Xn and the limit itself heavily depend on the sub- 
sequence Ui chosen to approach infinity. 

This should be compared with a more "regular" 
result in Theorem 3. Note that almost all the condi- 
tions of this theorem are satisfied in the diode model, 
except that here the distribution of In po is arithmetic 
(recall that the value In po = —00 is permissible), so 
it is the discreteness of the environment distribution 
that does not provide enough "mixing" and hence 
leads to such peculiar features of the asymptotics. 

9 Some Generalizations and 
Variations 

Most of the results discussed above in the simplest 
context of RWRE with nearest-neighbor jumps in an 
i.i.d. random environment, have been extended to some 
other cases. One natural generalization is to relax the 
i.i.d. assumption, e.g. by considering stationary er- 
godic environments (see details in Zeitouni (2004)). 
In this context, one relies on an ergodic theorem in- 
stead of the usual strong LLN. For instance, this way 
one readily obtains an extension of Solomon's cri- 
terion of transience vs. recurrence (see Theorem 1, 
Section 2). Other examples include an LLN (along 
with a formula for the asymptotic velocity, cf. The- 
orem 2, Section 3), a CLT and stable laws for the 
asymptotic distribution of X„ (cf. Theorem 3, Sec- 
tion 5), and Sinai's localization result for the recur- 
rent RWRE (cf. Theorem 4, Section 6). Usually, 
however, ergodic theorems cannot be applied directly 
(like, e.g., to X„, as the sequence Xn — Xn-i is not 
stationary). In this case, one rather uses the hitting 
times which possess the desired stationarity (cf. Sec- 
tions 3, 4). In some situations, in addition to station- 
arity one needs suitable mixing conditions in order to 
ensure enough decouphng (e.g., in Sinai's problem). 
The method of environment viewed from the parti- 
cle (see Section 7) is also suited very well to dealing 
with stationarity. 

In the remainder of this section, we describe some 
other generalizations including RWRE with bounded 
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jumps, RWRE where randomness is attached to bonds or equivalently 
rather than sites, and continuous-time (symmetric) 
RWRE driven by the randomized master equation. 



9.1 RWRE with Bounded Jumps 

The previous discussion was restricted to the case of 
RWRE with nearest-neighbor jumps. A natural ex- 
tension is RWRE with bounded iumps. Let L, R be 
fixed natural numbers, and suppose that from each 
site X G Z jumps are only possible to the sites x + i, 
i = —L, . . . ,R, with (random) probabilities 



R 



Px{i) > 0, ^ px{i) = 1. 



(35) 



We assume that the random vectors Px{-) determin- 
ing the environment are i.i.d. for different x G Z (al- 
though many results can be extended to the stationary 
ergodic case). 

The study of asymptotic properties of such a model 
is essentially more complex, as it involves products 
of certain random matrices and hence must use ex- 
tensively the theory of Lyapunov exponents (see de- 
tails and further references in Bremont (2004)). Lya- 
punov exponents, being natural analogs of logarithms 
of eigenvalues, characterize the asymptotic action of 
the product of random matrices along (random) prin- 
cipal directions, as described by Oseledec's multi- 
plicative ergodic theorem. In most situations, how- 
ever, the Lyapunov spectrum can only be accessed 
implicitly, which makes the analysis rather hard. 

To explain how random matrices arise here, let us 
first consider a particular case R = 1, L > 1. As- 
sume that Px{—L),px{l) > 5 > for all X G Z 
(ellipticity condition, cf. eqn (5)), and consider the 
hitting probabilities Un ■= Pni'^o < oo}, where 
To := min{i > : < 0} (cf. Section 2). By 
decomposing with respect to the first step, forn > 1 
we obtain the difference equation 



= Pn(l) Un+1 + ^^Pni-i) 



(36) 



i=0 



with the boundary conditions uq = ■ ■ ■ = U-l+i = 
1. Using that 1 = Pn{l) + Yli=oPn{-i), we can 
rewrite eqn (36) as 

L 

Pn(l) {Un - Un+l) = ^^Pni-i) {Un-i - Un) , 
1=1 



I Vn 



(37) 



where 



Ui 



i=l 

Ui+i and 



Kii) := + (38) 

PnW 

Recursion (37) can be written in a matrix form, Vn = 



MnVn-u where Vn := {vn 



Mn 



1 



■,Vn~L+l) 
■ ■ . hn{L) \ 





(39) 



V ... 1 / 

and by iterations we get (cf. eqn (10)) 

Vn = Mn--- MiFo, Vb = (1 - «l, 0, . . . , 0)^. 

Note that M„ depends only on the transition prob- 
ability vector and hence M„ ■ ■ ■ Mi is the prod- 
uct of i.i.d. random (non-negative) matrices. By 

Furstenberg-Kesten's theorem, the limiting behavior 
of such a product, as n ^ oo, is controlled by the 
largest Lyapunov exponent 



71 



lim n 



-1 



In \\M„ 



,Mi 



(40) 



(by Kingman's sub-additive ergodic theorem, the 
limit exists P-a.s. and is non-random). It follows 
that, Pq-sls., the RWRE X„ is transient if and only 
if 7i / 0, and moreover, X„ +oo (— oo) when 
7i < ( > 0), whereas limX„ = — oo, limX„ = 
-|-oo when 71 = 0. 

For orientation, note that if = p{i) are non- 
random constants, then 71 = In Ai, where Ai > is 
the largest eigenvalue of Mq, and so 71 < if and 
only if Ai < 1. The latter means that the character- 
istic polynomial ip{X) := dct(Afo — XI) satisfies the 
condition {-l)^(p{l) > 0. To evaluate det(Mo- 1), 
replace the first column by the sum of all columns 

and expand to get (p{l) = (-l)^"^(6i -I h 

Substituting expressions (38) it is easy to see that the 
above condition amounts to p{l) — > 
0, that is, the mean drift of the random walk is posi- 
tive and hence Xn +00 a.s. 

In the general case, L > 1, R > 1, similar con- 
siderations lead to the following matrices of order 
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d:=L + R-l (cf. eqn (39)) 

(an{R-l) ... a„(l) bn{l) ... hn{L)\ 

1 

1 

Mn = 

y 

where are given by (38) and 

/.X _ Pn(^)H ^Pn{R) 

"""^'^ Pn{R) 

Suppose that the ellipticity condition is satisfied in 
the form > (5 > 0, i 7^ 0, -L < i < i?, 

and let 71 > 72 > • • • > 7d be the (non-random) 
Lyapunov exponents of {M^}. The largest exponent 
71 is again given by eqn (40), while other exponents 
are determined recursively from the equalities 

71 + • • • + 7fe = lim n-^ In ||a'=(M„ • • • Mi)|| 

n-*oo 

(1 < < d). Here A denotes the external (anti- 
symmetric) product: x /\y = —yf\x (x,y G M*^), 
and a'^'M acts on the external product space A^'R'^, 
generated by the canonical basis {cj^ A • • ■ A Cj^. , 1 < 
ii < ■ ■ ■ < ik < d}, as follows: 

a''M{xi a • • • a Xfe) := M{xi) A • • • A M{xk). 

One can show that all exponents except 7^ are 
sign definite: 7^-1 > > 7i?+i. Moreover, it is 
the sign of 7^? that determines whether the RWRE is 
transient or recurrent, the dichotomy being the same 
as in the case R = 1 above (with 71 replaced by 
jr). Let us also mention that an LLN and CLT can 
be proved here (see Bremont (2004)). 

In conclusion, let us point out an alternative ap- 
proach due to Bolthausen and Goldsheid (2000) who 
studied a more general RWRE on a strip Z x { , 1 , . . . , 
m — I}. The link between these two models is given 
by the representation X„ = mYn + Zn, where m := 
max{L, R}, y„ e Z, Z„ e {0, . . . , m - 1}. Ran- 
dom matrices arising here are constructed indirectly 
using an auxiliary stationary sequence. Even though 
these matrices are non-independent, thanks to their 
positivity the criterion of transience can be given in 
terms of the sign of the largest Lyapunov exponent, 
which is usually much easier to deal with. An ad- 
ditional attractive feature of this approach is that the 
condition Px{R) > (P-a.s.), which was essential 
for the previous technique, can be replaced with a 
more natural condition F{px{R) > 0} > 0. 



9.2 Random Bonds RWRE 

Instead of having random probabilities of jumps at 
each site, one could assign random weights to bonds 
between the sites. For instance, the transition proba- 
biUties px = p{x, x + l,u>) can be defined by 

Cx,x+1 



Px 



Cx—l,x "I" Cx,x+1 



(41) 



where Cx,x+i > are i.i.d. random variables on the 
environment space 17. 

The difference between the two models may not 
seem very prominent, but the behavior of the walk in 
the modified model (41) appears to be quite different. 
Indeed, working as in Section 2 we note that 

Qx Cx—lx 



Px 



Px 



hence, exploiting formulas (11) and (41), we obtain. 



-a.s.. 



n-l 

E 



COI 



coi nE Cq/ — ^ 00, (42) 



since E Cq/ > 0. Therefore, /oo = 1, i.e. the random 
walk is recurrent (Po-a.s.). 

The method of environment viewed from the par- 
ticle can also be apphed here (see Sznitman (2004)). 
Similarly to Section 7, we define a new probabiUty 
measure Q = /{to) F using the density 

f{u) = Z"^(c_i,o(a;) + coiH), 

where Z = 2 E cqi is the normalizing constant (we 
assume that E cqi < 00). One can check that Q is in- 
variant with respect to the transition kernel eqn (41), 
and by similar arguments as in Section 7 we obtain 
that lim„_^oo Xn/n exists (Pg-a.s.) and is given by 

/ d{0, iu) q{du>) = Z-^ E [coi - c_i,o] = 0, 
Jci 

so the asymptotic velocity vanishes. 

Furthermore, under suitable technical conditions 
on the environment (e.g., cqi being bounded away 
from and 00, cf. eqn (5)), one can prove the fol- 
lowing CLT: 



lim Pq 

n— >oo 



<x}=^{x), (43) 



where = (Ecqi -Eco/)""^. Note that cr^ < 1 
(with a strict inequality if coi is not reduced to a con- 
stant), which indicates some slowdown in the spatial 
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spread of the random bonds RWRE, as compared to 
the ordinary symmetric random walk. 

Thus, there is a dramatic distinction between the 
random bonds RWRE, which is recurrent and dif- 
fusive, and the random sites RWRE, with a much 
more complex asymptotics including both transient 
and recurrent scenarios, slowdown effects and sub- 
diffusive behavior. This can be explained heuristi- 
cally by noting that the random bonds RWRE is re- 
versible, that is, m{x)p{x,y) = m{y)p{y,x) for 
all x,y G Z, with m{x) := Cx-\,x + Cx,x+\ (this 
property also easily extends to multidimensional ver- 
sions). Hence, it appears impossible to create ex- 
tended traps which would retain the particle for a 
very long time. Instead, the mechanism of the diffu- 
sive slowdown in a reversible case is associated with 
the natural variability of the environment resulting 
in the occasional occurrence of isolated "screening" 
bonds with an anomalously small weight Cx,x+i- 

Let us point out that the RWRE determined by eqn 
(41) can be interpreted in terms of the random con- 
ductivity model (see Hughes, 1996). Suppose that 
each random variable Cx^x+i attached to the bond 
(x, X + 1) has the meaning of the conductance of this 
bond (the reciprocal, c~;^^^, being its resistance). If 
a voltage drop V is applied across the system of 
successive bonds, say from to N, then the same 
current / flows in each of the conductors and by 
Ohm's law we have / = Cx^x+iVx^x+i, where Vx^x+i 
is the voltage drop across the corresponding bond. 
Hence 

TV JV 

which amounts to saying that the total resistance of 
the system of consecutive elements is given by the 
sum of the individual resistances. The effective con- 
ductivity of the finite system, cat, is defined as the 
average conductance per bond, so that 

1 ^ 

~ N ^ '-^x,x+l^ 
x=G 

and by the strong LLN, c^^ ^ Cq/ as AT — oo 
(P-a.s.). Therefore, the effective conductivity of the 
infinite system is given by c = (Ecqj^^) , and we 
note that c < E cqi if the random medium is non- 
degenerate. 

Returning to the random bonds RWRE, eqn (41), 
it is easy to see that a site j is recurrent if and only if 



the conductance Cj^oo between x and oo equals zero. 
Using again Ohm's law, we have (cf. eqn (42)) 

oo 
x=j 

and we recover the result about recurrence. 

9.3 Continuous-Time RWRE 

As in the discrete-time case, a random walk on Z 
with continuous time is a homogeneous Markov chain 
Xt,t G [0, oo), with state space Z and nearest neigh- 
bor (or at least bounded) jumps. The term "Markov" 
as usual refers to the "lack of memory" property, 
which amounts to saying that from the entire history 
of the process development up to a given time, only 
the current position of the walk is important for the 
future evolution while all other information is irrele- 
vant. 

Since there is no smallest time unit as in the dis- 
crete-time case, it is convenient to describe transi- 
tions of Xt in terms of transition rates characterizing 
the likelihood of various jumps during a very short 
time. More precisely, if Pxy{t) ■= P{Xt = y\Xo = 
x} are the transition probabilities over time t, then 
f or /i 

Pxy{h) = Cxyh + o{h) (x^y), 

Pxx{h) = l-hJ2(^xy + oih). (44) 

Equations for the functions Pxy{t) can then be de- 
rived by adapting the method of decomposition com- 
monly used for discrete-time Markov chains (cf. Sec- 
tion 2). Here it is more convenient to decompose 
with respect to the "last" step, i.e. by considering 
all possible transitions during a small increment of 
time at the end of the time interval [0, t -I- /i]. Using 
Markov property and eqn (44) we can write 

Cyx 

y^x 

+ POx{t)(l- hY^c^yj +o{h), 

which in the limit /i — > yields the master equation 
(or Chapman-Kolmogorov's forward equation) 

^POx{t) = V'lcj/xPOj/l*) - CxyPOx{t)}, 

y^x (45) 

Pox{0) = 6o{x), 
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where do{x) is the Kronecker symbol. 

Continuous-time RWRE are therefore naturally de- 
scribed via the randomized master equation, i.e. with 
random transition rates. The canonical example, orig- 
inally motivated by Dyson's study of the chain of 
harmonic oscillators with random couplings, is a sym- 
metric nearest-neighbor RWRE, where the random 
transition rates Cxy are non-zero only for y = x ± 1 
and satisfy the condition Cx,x+i = Cx+i,x, otherwise 
being i.i.d. (see Alexander et al. (1981)). In this case, 
the problem (45) can be formally solved using the 
Laplace transform, leading to the equations 



1-1 



s + G- + Gt = Q (XT^O), 
where G~, are defined as 

POx{s) -P0,x±l{s) 



Gx '■= Cx,x±l 



POx{s) 



(46) 
(47) 



(48) 



and pox{s) := J^Pox{t)e **dt. From eqs (47), 
(48) one obtains the recursion 



G 



1 1 

+ 



Cx,x±l S + G+j^i 

x = 0,±l,±2,... 



(49) 



The quantities are therefore expressed as infinite 
continued fractions depending on s and the random 
variables Cx.x±i, Cx.x±2, ■ ■ ■ The function poo(s) can 
then be found from eqn (46). 

In its generaUty, the problem is far too hard, and 
we shall only comment on how one can evaluate the 
annealed mean 

Epoo(s) = E(s + G+ + Go)-^ 

According to eqn (49), the random variables G^, Gq 
are determined by the same algebraic formula, but 
involve the rate coefficients from different sides of 
site X, and hence are i.i.d. Furthermore, eqn (49) 
implies that the random variables Gq, Gf have the 
same distribution and, moreover, Gf and cqi are in- 
dependent. Therefore, eqn (49) may be used as an 
integral equation for the unknown density function 
of Gq. It can be proved that the suitable solution 
exists and is unique, and although an explicit solu- 
tion is not available, one can obtain the asymptotics 
of small values of s, thereby rendering information 



about the behavior of poo (i) for large t. More specif- 
ically, one can show that if := (Kcq^) ^ > 
then 

Epoo(s) ~ (4c*s)-^/^ s^O, 



and so by a Tauberian theorem 



t — > oo. 



(50) 



Note that asymptotics (50) appears to be the same 
as for an ordinary symmetric random walk with con- 
stant transition rates Cx,x+i = Cx+i,x = c*, suggest- 
ing that the latter provides an "effective medium ap- 
proximation" (EMA) for the RWRE considered above. 

This is further confirmed by the asymptotic cal- 
culation of the annealed mean-square displacement, 
as t — ^ oo (Alexander et al. 1981). 
Moreover, Kawazu and Kesten (1984) proved that 
is asymptotically normal: 



lim P[ 



Xt 



< X 



(51) 



Therefore, if c^, > then the RWRE has the same 
diffusive behavior as the corresponding ordered sys- 
tem, with a well-defined diffusion constant D = c^. 

In the case where c* = (i.e., E Cq/ = oo), one 
may expect that the RWRE exhibits subdiffusive be- 
havior. For example, if the density function of the 
transition rates is modelled by 

fa{u) = (1 - a) u-"l{o<„<i} (0 < a < 1), 
then, as shown by Alexander et al. (1981), 

In fact, Kawazu and Kesten (1984) proved that in this 
case has a (non-Gaussian) Hmit distri- 

bution as t ^ oo. 

To conclude the discussion of the continuous-time 
case, let us point out that some useful information 
about recurrence of Xf can be obtained by consid- 
ering an imbedded (discrete-time) random walk X„, 
defined as the position of Xt after n jumps. Note 
that continuous-time Markov chains admit an alter- 
native description of their evolution in terms of so- 
journ times and the distribution of transitions at a 
jump. Namely, if the environment uj is fixed then 
the random sojourn time of Xt in each state x is ex- 
ponentially distributed with mean 1/cx, where Cx ■= 
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^y^x ^xy, while the distribution of transitions from 
X is given by the probabilities Pxy = Cxy/cx- 

For the symmetric nearest-neighbor RWRE con- 
sidered above, the transition probabilities of the imbed- 
ded random walk are given by 

Cx,x+1 

Px '■— Px,x+1 — ~ ) 

^x—l,x ~r Cx^x+1 

Qx '~ Px,x—1 ~ 1 Px: 

and we recognize here the transition law of a random 
walk in the random bonds environment considered 

in Section 9.2 (cf. eqn (41)). Recurrence and zero 
asymptotic velocity estabhshed there are consistent 
with the results discussed in the present section (e.g., 
note that the CLT for both X„, eqn (43), and Xt, 
eqn (51), does not involve any centering). Let us 
point out, however, that a "naive" discretization of 
time using the mean sojourn time appears to be in- 
correct, as this would lead to the scaling t = nSi 
with 5i := E(c_i^o + coi)~^, while from comparing 
the limit theorems in these two cases, one can con- 
clude that the true value of the effective discretization 
step is given by := (2c*)~^ = ^'KcqI. In fact, 
by the arithmetic-harmonic mean inequality we have 
5* > 6i, which is a manifestation of the RWRE's 
diffusive slowdown. 

10 RWRE in Higher Dimensions 

Multidimensional RWRE with nearest-neighbor 
jumps are defined in a similar fashion: from site x G 
the random walk can jump to one of the 2d ad- 
jacent sites X + e G Z*^ (such that |e| = 1), with 
probabilities (e) > 0, Y.\e\=iP'Ae) = 1, where 
the random vectors Px{-) are assumed to be i.i.d. for 
different x G Z^. As usual, we will also impose the 
condition of uniform ellipticity: 

Px{e)>d>0, |e| = 1, X G Z'', P-a.s.. (52) 

In contrast to the one-dimensional case, theory of 
RWRE in higher dimensions is far from maturity. 
Possible asymptotic behaviors of the RWRE for d > 
2 are not understood well enough, and many basic 
questions remain open. For instance, no definitive 
classification of the RWRE is available regarding tran- 
sience and recurrence. Similarly, LLN and CLT have 
been proved only for a limited number of specific 
models, while no general sharp results have been ob- 
tained. On a more positive note, there has been con- 
siderable progress in recent years in the so-called 



ballistic case, where powerful techniques have been 
developed (see Sznitman (2002, 2004) and Zeitouni 
(2003, 2004)). Unfortunately, not much is known for 
non-ballistic RWRE, apart from special cases of bal- 
anced RWRE ind > 2 (Lawler 1982), small isotropic 
perturbations of ordinary symmetric random walks 
in d > 3 (Bricmont and Kupiainen 1991), and some 
examples based on combining components of ordi- 
nary random walks and RWRE in d > 7 (Bolthausen 
et al. 2003). In particular, there are no examples of 
subdiffusive behavior in any dimension d > 2, and 
in fact it is largely beheved that a CLT is always true 
in any uniformly elliptic, i.i.d. random environment 
in dimensions d > 3, with somewhat less certainty 
about d = 2. A heuristic explanation for such a strik- 
ing difference with the case d= lis that due to a less 
restricted topology of space in higher dimensions, 
it is much harder to force the random walk to visit 
traps, and hence the slowdown is not so pronounced. 

In what follows, we give a brief account of some 
of the known results and methods in this fast devel- 
oping area (for further information and specific ref- 
erences, see an extensive review by Zeitouni (2004)). 

10.1 Zero-One Laws and LLNs 

A natural first step in a multidimensional context is 
to explore the behavior of the random walk X„ as 
projected on various one-dimensional straight lines. 
Let us fix a test unit vector £ G M'^, and consider the 
process := X„ • £. Then for the events A±( := 
{lim„_,.oo = ±00} one can show that 

Po{Ae U A_e) e {0, 1}. (53) 

That is to say, for each £ the probability that the ran- 
dom walk escapes to infinity in the direction £ is ei- 
ther or 1. 

Let us sketch the proof. We say that r is record 
time a \Zf\ > \Zl\ for all k < t, and regeneration 
time if in addition \Z^\ < \Z^\ for all n > r. Note 
that by the ellipticity condition (52), lim^^oo l-^nl = 
00 (Po-a.s.), hence there is an infinite sequence of 
record times = tq < ri < r2 < • • • If Po{A£ U 
A^fi) > 0, we can pick a subsequence of record 
times r/, each of which has a positive Pq -probability 
to be a regeneration time (because otherwise \Z^\ 
would persistently backtrack towards the origin and 
the event Ai U A^^ could not occur). Since the trials 
for different record times are independent, it follows 
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that a regeneration time r* occurs Pq-sls. Repeat- 
ing this argument, we conclude that there exists an 
infinite sequence of regeneration times r*, which im- 
plies that \Z^\ ^ oo (Po-a.s.), i.e., P{Ai U A_i) = 1. 

Regeneration structure introduced by the sequence 
{Tj* } plays a key role in further analysis of the RWRE 
and is particularly useful for proving an LLN and a 
CLT, due to the fact that pieces of the random walk 
between consecutive regeneration times (and frag- 
ments of the random environment involved thereby) 
are independent and identically distributed (at least 
starting from t^). In this vein, one can prove a "di- 
rectional" version of the LLN, stating that for each i 
there exist deterministic ve, V-g (possibly zero) such 
that 



lim =^ 

n— »oo n 



\ \a, + v-i \a_. , Po -a.s. (54) 



Note that if PoiAe) G {0, 1}, eqn (54) in con- 
junction with eqn (53) would readily imply 



lim — = vi, Po-a.s. 

n— »oo n 



(55) 



Moreover, if Po^Ag) G {0, 1} for any i, then there 
exists a deterministic v (possibly zero) such that 



iim — 



V, Po -a.s. 



(56) 



Therefore, it is natural to ask if a zero-one law (53) 
can be enhanced to that for the individual probabili- 
ties Poi^Ai). It is known that the answer is affirma- 
tive for i.i.d. environments in d = 2, where indeed 
P{Ae) G {0, 1} for any £, with counter-examples in 
certain stationary ergodic (but not uniformly elliptic) 
environments. However, in the case d > 3 this is an 
open problem. 

10.2 Kalikow's Condition and Sznitman's 
Condition (T') 

An RWRE is called ballistic (ballistic in direction £) 
if V ^ (vi ^ 0), see eqs (55), (56). In this section, 
we describe conditions on the random environment 
which ensure that the RWRE is ballistic. 

Let U he a. connected strict subset of Z*^ contain- 
ing the origin. For x e U, denote by 

Tu 
n=0 



the quenched mean number of visits to x prior to the 
exit time Tu := min{n > : ^ U}. Con- 
sider an auxiliary Markov chain Xn, which starts 
from 0, makes nearest-neighbor jumps while in U, 
with (non-random) probabilities 



Px{e) 



E[g{x,uj)p^{e)] 
E[g{x,io)] ■ 



X eU, 



(57) 



and is absorbed as soon as it first leaves U. Note that 
the expectations in eqn (57) are finite; indeed, if ax is 
the probability to return to x before leaving U, then, 
by the Markov property, the mean number of returns 
is given by 



^A;a^(l - a. 



a.7 



k=l 



< oo, 



since, due to ellipticity, a.^ < I. 

An important property, highlighting the usefulness 
of Xn, is that if Xn leaves U with probability 1, then 
the same is true for the original RWRE X„ (under 
the annealed law Po), and moreover, the exit points 
Xf^ and Xt^ have the same distribution laws. 

Let i G M'', 1^1 = 1. One says that Kalikow's con- 
dition with respect to I holds if the local drift of Xn 
in the direction I is uniformly bounded away from 
zero: 

inf inf ^ (e • ^) > 0. (58) 

A sufficient condition for (58) is, for example, that 
for some k> 

E[{d{0,uj) ■£)+] > KE[{d{0,uj) ■£)-], 

where d(0, uj) = Xi and u± := maxjztit, 0}. 

A natural implication of Kalikow's condition (58) 
is that PoiA^) = 1 and > (see eqn (55)). More- 
over, noting that eqn (58) also holds for all i' in a 
vicinity of £ and applying the above result with d 
non-coUinear vectors from that vicinity, we conclude 
that under Kalikow's condition there exists a deter- 
ministic V such that Xn/n ^ u as n ^ oo 
(Po-a.s.). Furthermore, it can be proved that {Xn — 
nv) I ^Jn converges in law to a Gaussian distribution 
(see Sznitman (2004)). 

It is not hard to check that in dimension d = \ 
Kalikow's condition is equivalent to f / and there- 
fore characterizes completely all ballistic walks. For 
d > 2, the situation is less clear; for instance, it is 
not known if there exist RWRE with P{Ai) > and 
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ve = (of course, such RWRE cannot satisfy Ka- 
likow's condition). 

Sznitman (2004) has proposed a more compUcated 
transience condition (T') involving certain regenera- 
tion times T* similar to those described in Section 
10.1. An RWRE is said to satisfy Sznitman' s condi- 
tion (T') relative to direction I if P^i^Ai) = 1 and 
for some c > and all < 7 < 1 

-Egexplc sup l-^nf' ) < oo. (59) 

^ n<T* ^ 

This condition provides a powerful control over rf 
for d > 2 and in particular ensures that t\ has finite 
moments of any order. This is in sharp contrast with 
the one-dimensional case, and should be viewed as a 
reflection of much weaker traps in dimensions d > 
2. Condition (59) can also be reformulated in terms 
of the exit distribution of the RWRE from infinite 
thick slabs "orthonormal" to directions £' sufficiently 
close to £. As it stands, the latter reformulation is 
difficult to check, but Sznitman (2004) has developed 
a remarkable "effective" criterion reducing the job 
to a similar condition infinite boxes, which is much 
more tractable and can be checked in a number of 
cases. 

In fact, condition (T') follows from Kalikow's con- 
dition, but not the other way around. In the one- 
dimensional case, condition (T') (applied to ^ = 1 
and i = —1) proves to be equivalent to the transient 
behavior of the RWRE, which, as we have seen in 
Theorem 2 (Section 3), may happen with v = 0, i.e. 
in a non-ballistic scenario. The situation in d > 2 
is quite different, as condition (T') implies that the 
RWRE is ballistic in the direction £ (with > 0) and 
satisfies a CLT (under Pq). It is not known whether 
the ballistic behavior for d > 2 is completely char- 
acterized by condition (T'), although this is expected 
to be true. 

10.3 Balanced RWRE 

In this section we discuss a particular case of non- 
ballistic RWRE, for which LLN and CLT can be proved 
Following Lawler (1982), we say that an RWRE is 

balanced if ^^^(e) = — e) for all x G Z"^, |e| = 1 
(P-a.s.). In this case, the local drift vanishes, d{x, u) = 
0, hence the coordinate processes (i = 1,. . . ,d) 
are martingales with respect to the natural filtration 
= (^{Xq, . . . , Xn}- The quenched co variance 
matrix of the increments AX^ := X!.,^ — X^ {i = 



1 , . . . , d) is given by 

E'^ [AXi AXi I Tr^] = •^SijPxM). (60) 

Since the right-hand side of eqn (60) is uniformly 
bounded, it follows that Xn/n (Po-a.s.). Fur- 
ther, it can be proved that there exist deterministic 
positive constants ai , . . . , such that for z = 1 , . . . , d 

^ n— 1 

lim - Vpxfc(ei) = Po-a.s. (61) 

fc=0 

Once this is proved, a multidimensional CLT for mar- 
tingale differences yields that X^./ ^Jn converges in 
law to a Gaussian distribution with zero mean and 
the CO variances hij = bi^ai. 

The proof of (61) employs the method of environ- 
ment viewed from the particle (cf. Section 7). Namely, 
define a Markov chain a;„ := Q-^^uj with the transi- 
tion kernel 

d 

i=l 

+ Po{-ei) 5f)-i^{duj')] 

(cf. eqn (27)). The next step is to find a probabil- 
ity measure Q on O invariant under T and absolutely 
continuous with respect to P. Unlike the one-dimensional 
case, however, an explicit form of Q is not available, 
and Q is constructed indirectly as the limit of invari- 
ant measures of certain periodic modifications of the 
RWRE. Birkhoff 's ergodic theorem then yields, Po- 
a.s., 

^ n— 1 ^ n— 1 

n ^-^ n ^-^ 

k=o k=o 

/ po{ei,u;i)Q{dLo) > S, 
Jn 

by the elUpticity condition (52), and eqn (61) fol- 
lows. 

With regard to transience, balanced RWRE admit 
a complete and simple classification. Namely, it has 
been proved (see Zeitouni (2004)) that any balanced 
RWRE is transient for d > 3 and recurrent for d = 2 
(Po-a.s.). It is interesting to note, however, that these 
answers may be false for certain balanced random 
walks in a fixed environment (P-probability of such 
environments being zero, of course). Indeed, exam- 
ples can be constructed of balanced random walks in 
and in Z** with d > 3, which are transient and 
recurrent, respectively (Zeitouni 2004). 
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10.4 RWRE Based on Modification of 
Ordinary Random Walks 

A number of partial results are known for RWRE 
constructed on the basis of ordinary random walks 
via certain randomization of the environment. A nat- 
ural model is obtained by a small perturbation of a 
simple symmetric random walk. To be more precise, 
suppose that: (a) \px{e) — ^\ < £ for all x G Z*^ 
and any |e| = 1, where e > is small enough; 
(b) Epx{e) = (c) vectors Px{-) are i.i.d. for dif- 
ferent X £ 'Z'^, and (d) the distribution of the vec- 
tor Px{ ) is isotropic, i.e. invariant with respect to 
permutations of its coordinates. Then for d > 3 
Bricmont and Kupiainen (1991) have proved an LLN 
(with zero asymptotic velocity) and a quenched CLT 
(with non-degenerate covariance matrix). The proof 
is based on the renormahzation group method, which 
involves decimation in time combined with a suit- 
able spatial-temporal scaling. This transformation 
replaces an RWRE by another RWRE with weaker 
randomness, and it can be shown that iterations con- 
verge to a Gaussian fixed point. 

Another class of examples are also built using small 
perturbations of simple symmetric random walks, but 
are anisotropic and exhibit ballistic behavior, pro- 
viding that the annealed local drift in some direction 
is strong enough (see Sznitman (2004)). More pre- 
cisely, suppose that d > 3 and t] G (0, 1). Then there 
exists £0 = £o{d, '?) > such that if \px{e) ~ ^ I < ^ 
(x G Z*^, |e| = 1) with < e < eo, and for some eo 
one has E [d{x, oj) ■ eo] > e^-^-'> {d = 3) or > e^-'? 
{d > 4), then Sznitman's condition (T') is satisfied 
with respect to eo and therefore the RWRE is ballis- 
tic in the direction eo (cf. Section 10.2). 

Examples of a different type are constructed in di- 
mensions d > 6 by letting the first di > 5 coordi- 
nates of the RWRE Xn behave according to an ordi- 
nary random walk, while the remaining d2 = d — di 
coordinates are exposed to a random environment 
(see Bolthausen et al. (2003)). One can show that 
there exists a deterministic v (possibly zero) such 
that Xn/n V (Po-a.s.). Moreover, if c/i > 13 
then {Xn — nv) / ^Jn satisfies both quenched and an- 
nealed CLT. Incidentally, such models can be used 
to demonstrate the surprising features of the multidi- 
mensional RWRE. For instance, for d > 7 one can 
construct an RWRE X„ such that the annealed local 
drift does not vanish, E d(x, uS) / 0, but the asymp- 
totic velocity is zero, X^jn (Po-a.s.), and fur- 



thermore, if d > 15 then in this example Xn/^/n 
satisfies a quenched CLT. (In fact, one can construct 
such RWRE as small perturbations of a simple sym- 
metric walk.) On the other hand, there exist exam- 
ples (in high enough dimensions) where the walk is 
ballistic with a velocity which has an opposite di- 
rection to the aimealed drift 'Kd{x,uj) ^ 0. These 
striking examples provide "experimental" evidence 
of many unusual properties of the multidimensional 
RWRE, which, no doubt, will be discovered in the 
years to come. 

See also: Averaging Methods; Growth Processes in 
Random Matrix Theory; Lagrangian Dispersion 
(Passive Scalar); Random Dynamical Systems; Ran- 
dom Matrix Theory in Physics; Stochastic Differen- 
tial Equations; Stochastic Loewner Evolutions. 
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